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In this review article, we survey the relatively new theory of orbital magnetization in 
solids — often referred to as the "modern theory of orbital magnetization" — and its appli- 
cations. Surprisingly, while the calculation of the orbital magnetization in finite systems 
such as atoms and molecules is straight forward, in extended systems or solids it has 
long eluded calculations owing to the fact that the position operator is ill-defined in such 
a context. Approaches that overcome this problem were first developed in 2005 and in 
the first part of this review we present the main ideas reaching from a Wannier function 
approach to semi-classical and finite-temperature formalisms. In the second part, we de- 
scribe practical aspects of calculating the orbital magnetization, such as taking fc-space 
derivatives, a formalism for pscudopotentials, a single fc-point derivation, a Wannier in- 
terpolation scheme, and DFT specific aspects. We then show results of recent calculations 
on Fe, Co, and Ni. In the last part of this review, we focus on direct applications of the 
orbital magnetization. In particular, we will review how properties such as the nuclear 
magnetic resonance shielding tensor and the electron paramagnetic resonance g-tensor 
can elegantly be calculated in terms of a derivative of the orbital magnetization. 

Keywords: Magnetism; orbital magnetic moment and magnetization; solids; periodic and 
extended systems; Berry phase; Wannier functions; Bloch functions; semi-classical; finite 
temperature; pseudopotentials; Chern invariant; nuclear magnetic resonance; electron 
paramagnetic resonance; first-principles; density functional theory. 

1. Introduction 

Magnetism is one of the oldest phenomena known to mankind. In the forms of lode- 
stone and magnetite it is first mentioned as early as 4000 BC in Chinese writings. 1 
The name is often related to the Greek province of Magnesia where the mineral 
was first mined, or to the shepherd Magnes, "the nails of whose shoes and the 
tip of whose staff stuck fast in a magnetick field while he pastured his flocks. '00 
However, it was not until the 19 th century that scientists such as Oersted, Ampere, 
Gauss, Faraday, and Maxwell took much of the mysticism out of magnetism and 
put its theory on a stronger, classical scientific foundation. Magnetism's quantum- 
mechanical origin was then developed during the first half of the 20 th century. Today, 
many standard textbooks on the subject exist, courses in Electricity and Magnetism 
are required for many undergraduate and graduate students, and magnetism is om- 
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nipresent in today's life through a vast amount of technical applications and gadgets 
that rely on it. As such, it is rather surprising that, until just recently, no complete 
fundamental theory for magnetism in solids existed. 

Nowadays, it is well understood that magnetism in atoms and molecules orig- 
inates from two distinct sources related to quantum- mechanics: the spin magnetic 
moment and the orbital magnetic moment of the electrons^] Rigorous theories for 
both parts exist at several levels, and thus, the calculation of magnetic moments 
of atoms and molecules has become a standard procedure. The situation in solids 
is similar; the quantity of interest here is the magnetization, i.e. the magnetic mo- 
ment per volume. In parallel to atoms and molecules, the magnetization in solids 
has two distinct contributions: the spin magnetization and the orbital magnetiza- 
tion. However, while the first can be computed easily, the latter has eluded exact 
calculation due to the lack of an appropriate theory. It was not until 2005 that 
the corresponding theory was developed. Note that already for several decades the 
two contributions can be determined separately through a combination of magneto- 
mechanical measurements Pin the present review article we focus on this relatively 
new modern theory of orbital magnetization in solids, collecting the fundamental 
derivations, its applications, and first results. 

The focus of this review article is entirely on the orbital contribution to mag- 
netism in solids — as the title emphasizes. For the purpose of this review article we 
may refer to solids also as infinite, periodic, or extended systems, or simply as crys- 
tals. The calculation of the orbital contributions in finite systems, i.e. atoms and 
molecules, is well established and is not the subject of this review; it can readily be 
found in standard textbooks. 

This review article is organized in the following way: In Section [2j we explain 
the importance of the orbital magnetization and its relation to several other proper- 
ties of current interest. Section [3] contains the theoretical framework for the theory 
of orbital magnetization. First, in Section [3T| we focus on the difference between 
finite systems and solids and explain the mathematical difficulties that arise in the 
context of periodic systems. Then, we show several different derivations of the main 
theory, from a Wannier function approach in Sections |3 . 2 1 and [3 . 3| to a semi-classical 



(Section 3.4) and finite-temperature (Section 3.5 I formalism, which can be extended 
to interacting systems. Practical aspects of calculating the orbital magnetization in 
solids are collected in Section [4] In particular, we review options for taking the 
required fc-spacc derivatives in Section 4.1 and how to perform calculations using 
pseudopotentials in Section 4.2 We further review a single fc-point formula (Sec- 



tion |4^3|) , a Wannier interpolation scheme (Section 4.4), and DFT specific aspects 
(Section 4.5). In Section [5] we summarize first results for calculations of the or- 
bital magnetization in Fe, Co, and Ni. Finally, we review several applications of the 
orbital magnetization in Section[6j In particular, Sections 6.1 and 6.2 focus on show- 



a Here, we are neglecting the magnetic moment originating from the nuclear spins, which are usually 
smaller by several orders of magnitude. 
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ing how theories of nuclear magnetic resonance (NMR) and electron paramagnetic 
resonance (EPR) can elegantly be expressed as derivatives of the orbital magneti- 
zation. Other important derivatives are reviewed in Section [6.3| We conclude with 
an outlook in Section 



2. Importance of the Orbital Magnetization 

As pointed out in the introduction, magnetism in solids originates from two distinct 
sources, i.e. spin magnetic moments and orbital magnetic moments. As such, the 
orbital magnetization is a crucial part of the overall magnetization. The reason 
why it has not attracted more attention in the past — and also the reason why 
the condensed-matter/materials-science community "survived" without a theory for 
it — is that in many common materials of everyday interest, the orbital contribution 
is small compared to the spin contribution. In Fe, Co, and Ni for example, it is 
only of the order of a few percent of the total magnetizationP not surprising, its 
effects have been neglected in the past. However, that is not to say that the orbital 
magnetization is small in all materials and in other systems it can have a more 
important effect and even cancel the spin magnetization.^HH Overall, the orbital 
magnetization is important in-and-of-itself for a complete description of magnetism 
in solids. 

Furthermore, a wealth of applications are directly related to the orbital magne- 
tization. For example, a theory of solid-state NMR can conveniently be defined in 



terms of a derivative of the orbital magnetization^ discussed in detail in Section 6.1 
removing difficulties related to gauge-origin choices of previous methods. In a simi- 
lar sense, the EPR g-tensor follows from a derivative of the orbital magnetization^ 
as described in Section [672] Other derivatives relate to properties such as the mag- 
netic susceptibility, the orbital magnetoelectric coupling and responsej^tni the spin 
Hall conductivity,^ and the identification of non-abelian quantum Hall states.^ 

The modern theory of orbital magnetization is further important because of its 
close connection to the modern theory of electric polarization in solids J 1 H 15 l In fact, 
there is an underlying connection between these two theories through the physics of 
the Berry phased which has been discovered to thread through many different fields 
of physics. This connection between Berry's phase, the electric polarization, and the 
orbital magnetization is thoroughly discussed in two other review articles EZEH 

Orbital magnetization permeates all dimensions — it plays an important role in 
zero-dimensionalf^mS] one-dimensionalp3H22l two-dimensionalJ^^H an d obviously 

three-dimensional systems HEEMU] The orbital magnetization is of interest in 
adatoms ji£IE2l clusters and quantum dotspH^ID and in metal- and nano- wires p3HHI \i 
plays a role in graphene^^ and the two-dimensional electron gasP^SU ft further is 
significant for certain ferromagnets^ and antiferromagnetsf^ in some ferromagnets 
it can even compensate the spin magnetization, resulting in zero net magnetiza- 
tionPH The orbital magnetization is of importance in semiconductors,^ metallic 
magnets,^ itinerant magnets,^ and superconductorsE2E3 The orbital magnetiza- 
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tion can play an important part in certain half- met alsp^ Mott insulators,^ and 
ordered alloysP^It influences the thermoelectric, anomalous, and topological trans- 
port jHllMllII] the magnetic susceptibility,^ as well as the anisotropic magnetoresis- 
tanceP^ Finally, it also shows interesting effects at interfaces,^ in the quantum Hall 
fluid in localized electrons,^ and in connection with the orbital polarization of 
itinerant magnetsP^ 



3. Theory of Orbital Magnetization in Solids 

The theory of orbital magnetization was developed by two research groups indepen- 
dently in 2005. Although their approaches are completely different, they complement 
each other and the final expressions derived within each approach — which we have 
indicated with boxes for easier reference — all come to agreement. In the following 
sections we review the theory presented in the corresponding original manuscripts. 
Note that we are consistently using Gaussian units throughout, which, at times, 
deviates from the original manuscripts. However, before showing details about the 
derivation of the theory of orbital magnetization in solids, it is instructive to take 
a moment and consider the original problem. 



3.1. What's the problem? 

As pointed out in the introduction, the calculation of orbital contributions to mag- 
netism in atoms and molecules is rather straightforward. But, what makes its cal- 
culation so difficult in solids? 

If we consider a finite system, the orbital contribution to the magnetic moment 
can be expressed in terms of basic quantities. In a simplified one-particle picture, 
defining the velocity operator v as the commutator of the position operator r with 
the system's Hamiltonian H, 

v = - l ~[r,H], (1) 
the total orbital moment in Gaussian units becomes 

m orb = -— ^/ra (">Pn\* X V|^n> • (2) 
n 

Here, \tp n ) & r e the eigenstates of H, f n is the occupation number, and — e and c are 
the electronic charge and the speed of light, respectively. This is in parallel to the 
classically defined moment caused by a current density J(r) in Gaussian units,^ 

m orb = i J d 3 r r x J = i J d 3 r r x pw , (3) 

where in the last step we have explicitly included the charge density p(r). The 
magnetization M or b can then be defined as the magnetic moment per unit volume 



Theory of Orbital Magnetization in Solids 5 



V as0 



M orb = ^ = E/« ^«l r x v ^> • W 



Equations |2]) and Q are valid for all finite systems, even when very large. 

The essential difficulty for truly periodic systems becomes apparent when taking 
the thermodynamic limit of Eq. Q. In such systems, the eigenstates \ip n ) become 
Bloch states \ipnk) with Bloch wave vector k, in-between which the position operator 
is ill-defined since Bloch states are extended. It is tempting to try to solve the 
problem by simply changing to a more localized basis, but this alone does not solve 



the problem, as we will show in Section 3.2 Also, linear-response methods exist for 



periodic systems that circumvent the problem in other ways, but they only allow 
the calculation of magnetization changes, not of the magnetization itself.^HSU 

One might also hope to solve the problem from the point of view of the local 
bulk current density J(r) and its relation to the magnetization. However, while the 
eigenstates \i/j n k) can be used to calculate J(r), it is often wrongly assumed that 
this uniquely defines the magnetization. As Hirst has emphasized,^ the knowledge 
of J(r) — even in principle — is insufficient to calculate the macroscopic orbital mag- 
netization M or b- This can easily be seen from the following argument. Consistent 
with Maxwell's equations we can define the orbital magnetization density Al rb(r) 
asP 

cV x M orb (r) = J(r) . (5) 

However, A4 rb(r) can simply be replaced by another .A/t^^r), i.e. 

M orb (r) M orb (r) = M orb (r) + M° olb + V£(r) , (6) 

which has in addition a constant shift M-^ib an d even the gradient of an arbitrary 
function £(r), but nevertheless produces the same local bulk current density J(r)PSl 
Hence, the macroscopic orbital magnetization M orb cannot be uniquely defined as 
the cell average of A / t orb (r). In this respect, the theory of orbital magnetization is 
similar to the theory of electric polarization, where knowledge of the density p(r) — 
even in principle — is insufficient to determine the electric polarization PlHILLMH] 

Since the calculation of the orbital magnetization in solids proves difficult, a 
simple approximation is often used, referred to as muffin-tin approximation. Non- 
overlapping muffin-tin spheres can be centered around the atoms in a solid. Within 
these spheres, which are finite systems, the moment can then be calculated according 
to Eq. pi). Often, the orbital magnetization indeed originates from regions near 
the atom cores, making this approximation good when the magnetization is mostly 
confined to the muffin-tin spheres and the interstitial contribution can be neglected. 



This muffin-tin approximation has for example been used in Refs. 59 60 . Now 



b The definition of a magnetization, i.e. a magnetic moment per unit volume, is formally only 
useful for periodic systems, but we use it here for finite systems representing a fragment of a 
larger, crystalline system. 
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having the full theory of orbital magnetization, we can assess the accuracy of this 
approximation quantitatively, as done for Fe, Co, and Ni in Section [5j 

3.2. Wannier function derivation of the orbital magnetization 

For this first derivation of the orbital magnetization we limit ourselves to insulating 
solids, described by a one-particle Hamiltonian with broken time-reversal symme- 
try. We further restrict ourselves to spinless electrons, one (completely occupied) 
bulk band, and vanishing Chern invariant These restrictions are non-essential 
and can easily be removed; a corresponding multi-band case, metallic systems, and 
Chern insulators are discussed in Section [3.3| By assuming either a vanishing macro- 
scopic field or an integral number of flux quanta per unit cell, the Bloch wave vector 
k is guaranteed to be a good quantum number. Breaking of time-reversal symmetry 
in such situations can occur through spin-orbit coupling to underlying ordered local 
moments )S2i2Sl We only focus on the main concepts of the Wannier function deriva- 



tive here and refer the reader to the original manuscript in Ref. 66 for further 
details. 

We start by considering a finite sample, which we then make larger and larger 
until it becomes a periodic solid. The orbital magnetization of the finite sample can 



be calculated using Eq. Q. However, as discussed in Section 3.1 the thermodynamic 



limit of this expression is not well-defined. In parallel to the situation of the electric 
polarization, it is tempting to switch from the delocalized Bloch-like eigenstates 
\ip n ) in Eq. Q, which cause the problem, to well-localized orthonormal orbitals 
\4> n ), simply by using the invariance of the trace. In the thermodynamic limit, the 
\4> n ) deep inside the sample can then be associated with bulk Wannier functions 
|R), where R is a lattice vector. Since Wannier functions are exponentially localized 
in insulatorspH^H matrix elements of the position operator r are now indeed well 
defined. It is at this point that this derivation relies on the restrictions described 
above, as Wannier functions in metals or Chern insulators are not exponentially 
localized 

Surprisingly, simply evaluating Eq. Q for a solid in-between the bulk Wannier 
function |R) in the home unit cell, i.e. |0), does not result in the total orbital 
magnetization!^ Instead, this leads to a contribution that corresponds to the locat- 
or self-circulation of a Wannier function 

M^ = --^-<0|rxv|0), (7) 
2cVb 

i.e. the circulation resulting from a Wannier function rotating around its center. 
The volume Vq corresponds to the unit-cell volume. Note that the sum symbol 
from Eq. Q disappears because we have limited ourselves to only one occupied 
bulk band, resulting in only one occupied Wannier function. Interestingly, the same 
approach in the case of the electric polarization P, i.e. replacing \ip n ) in the analog 
of Eq. Q with bulk Wannier functions (ip n \r\ip n ) — » (0|r|0), solved the problem 
completelyEEllim 
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Fig. 1. a) lOx 10 sample of the Haldane modelESlThe currents 
Wannier function are plotted usini 



e(<j) n \v\<j> n ) associated with each 
arrows at the corresponding Wannier centers r n = (<f> n \r\4>n) ■ 

to the orbital magnetization. 

b) Horizontal slice of the same materials as in Fig.JI^,, but infinitely long in the vertical y-direction. 
(Reprinted with permission from Ref. 166 ; © 2005 American Physical Society). 



The itinerant current gives rise to an additional contribution M'SL 



The missing, second contribution to the magnetization arises from itinerant sur- 
face contributions Mj~ b , which — even in periodic systems — have a remnant, fully 
contained in the bulk band structure. To find this missing contribution, we need to 
find the surface current caused by the itinerant circulation of the Wannier functions. 
Note that bulk Wannier functions of the sample, as they correspond to fully occu- 
pied bulk bands, cannot carry any current and thus cannot produce an itinerant 
contribution. Surface Wannier functions, on the other hand, behave differently; they 
are not eigenstates of the Hamiltonian and thus evolve over time. It is important 
to note that this surface contribution is different from the simple textbook picture 
where an array of current loops cancel inside a domain and only the surface cur- 
rent remains. This itinerant contribution originates from the fact that the surface 
Wannier functions or current loops themselves move around the surface! 

To show this fact, the authors of the original article in Ref. [66] used a sim- 
ple, two-dimensional tight-binding model. 61 A finite system consisting of 10x10 
unit cells of this hexagonal system is shown in Fig. [IJl This figure also shows re- 
sults for the calculation of the itinerant current — ev„ = — e(<fi n \v\cj) n ) (indicated 
by arrows) associated with each Wannier function, plotted at the Wannier centers 
?n = (4>n\r\4>n) ■ The current dies-off exponentially into the bulk, as bulk Wannier 
functions do not carry any current. The surface current obviously gives an additional 
contribution to the magnetization according to 

M - b = ~2cV^ f " X( ^ |v|0 " ) - (8) 

n 

Below, we show how this itinerant surface contribution can be expressed as a 
bulk property. To this end, we focus on the right edge of the sample displayed in 
Fig. [T}?. We use labels s and s' to denote Wannier functions to the left and right of 
the boundary at Xq, respectively. The y-component of the current on the right edge 
of the sample then becomes 

s' s' s s' 
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where we have used Eq. (JlJ and inserted a complete set of Wannier func- 
tions in-between r and H. In the last step we used the notation V/jj) — 
(2/H) Im ((f>i\r\<fij}(<f)j\H\(j)i}. Note that only Wannier functions that are within the 
vertical extension Al of the strip contribute to the primed sum. In principle, the 
first sum should run over all Wannier functions n, but the ones where n is within 
the subset of s' do not contribute to the overall current. Moving the boundary xq 
deep into the bulk and using translational symmetry, the current in Eq. ([9| becomes 

Iy = ~Al E E v (0,w-B.), y = ~^Y RxV (°' R )-y ' ( 10 ) 

Rx<x R' m >xo R 

where we have further used the fact that there arc exactly (R' x — R x )AI/Vq terms in 
the sum with a given value of R' — R if R' x — R x > and zero otherwise. We divide 
by 2 in order to extend the sum over all R. The z-component of the corresponding 
magnetic moment then finally becomes 

Mor C b,, = ~2c^y Im E (^(RMO)(0|H|R) - ^(R|r,|0)<0|ff|R)) , (11) 



■o 

or more generally 



R 



M - b = -4cF E Rxv <°,*>< ( 12 ) 

R 

clearly expressing also the itinerant contribution M^ b as a bulk property. 

The final step is to transform the expressions for and M Q ^ b back to the 

Bloch representation. With the definition of Wannier functions in terms of the cell- 
periodic part of the Bloch functions \u k ) = e~ l^k), i-e. 

I R ) = 7^ / d3k e 2k ' (r - R) |u k ) , (13) 



(14) 



(2tt) 3 d 

the main result becomes, after some straightforward but tedious algebra 



M orb = M££ + M D ? b = — ImJ — (cW| x (H k + E w ) \d k u k ) 



Here, we have used the shortcuts &k = <9/<9k and H k = e~ r H e ik r ; the E k denote 
the corresponding eigenvalues. This expression is reminiscent of a result found earlier 
for the special case of the Hofstadtcr model.^1 



Equation ( 14 ) is the main result of this section. It gives the orbital magnetiza- 
tion of a solid in terms of a Brillouin-zone integral over well- understood quantities. 
Even in the final expression, the two contributions corresponding to the internal 
circulation of bulk Wannier functions (containing H k ) and the net currents carried 
by Wannier functions near the surface (containing E k ), are clearly distinguishable. 
Note that the surface contribution can be written as 

M - = -2s/(Si E - n -- < 15 > 
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where Jlk = Vk x Ak is the Berry curvature and Ak = i(itk| Vk|ltk) the Berry 
connection — revealing the underlying relationship with the theory of the Berry 
phaseP^ffllD it [ s interesting to note that both contributions to the orbital mag- 
netization in Eq. ( 14 1 are individually gauge invariant^ and thus are, in principle, 
separately observable f^2l However, this is only true for the simplified case of a single 
occupied bulk band. 

Although maximally localized Wannier functions were used to derive Eq. ( 14 1, in 
the end the orbital magnetization can be written in terms of Bloch states — already 
suggesting that this expression might also be true for systems where maximally 
localized Wannier functions cannot be constructed, such as metals and Chern insu- 
lators. 70 In fact, in the next sections we will see that this is indeed the case. 



3.3. Multi-band derivation of the orbital magnetization 
and extensions to metals and Chern insulators 

The derivation of the orbital magnetization in the previous section was limited in 
that it only applied to insulators with vanishing Chern invariant and one occupied 
bulk band. Not long after the original theory appeared in Rcf 
were removed in Rcf. 



(>(» 



73 



these limitations 
and refer 



73 



We review here the main results of Ref. 
the reader to the original manuscript for further details. 

The derivative of a theory of orbital magnetization for multi-band systems fol- 
lows, in spirit, the outline of the previous section. We now assume that there 
are several occupied bulk bands in the solid, indicated by the band index n in 
|"0nk) = e lk r |w,ik), with occupation /„k- This leads to the existence of several Wan- 
nier functions in each unit cell labeled by |nR). It is easy to see that, following 
the same arguments from Section |3.2[ one arrives at the same equations, with the 
only difference being that the Wannier functions get an index and sums over bands 
occur. In summary, Eqs. ^ and ( |12[ ) turn into 

(nOlr x vlnO) , 



M L K = - 

■ LVJ -orb 



■ LVJ -orb 



2cV v 

— y 

nmtt 



R X V (mO,nR) 



(16) 
(17) 



which, after transformation back to Bloch states, becomes 



- LV - L orb 



M 



ic 



2hc ^ 



d 3 k 



/nk (<9 k Unk| X (# k + E n u) |<9k"U„ k ) 



(18) 



The main difference between Eq. (18) and Eq. (14) — besides the band indices 



and sum over bands — is the fact that in the single-band case both contributions 



c Here and in the following, we use the term gauge invariance to refer to the fact that a certain 
property is independent of the choice of phase factors of the Bloch states. All physical observablcs 
have to be gauge invariant in this sense. 
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M-orb and M^ b by themselves are gauge invariant, while in the multi-band case 
only their sum is. The corresponding proof goes beyond the scope of this review, 



but is given in much detail in Ref. 73 . In the same reference, it is also shown that 



terms in Eq. ( 18 ) can be regrouped such that the overall magnetization is again 



written as a sum of two contributions, whereas each contribution itself is now gauge 
invariant. This relationship is further investigated in Ref. [74] , where it is shown — 
through an elegant relation to the /-sum rule — that the two contributions can be 
determined separately through a combination of gyromagnetic and magneto-optical 
experiments. 



Similar to the single-band case below Eq. (15), one can define a multi-band 
Berry curvature as 



O nk = -i(d k u nk \ x |<9kti„ k ) , 
which in turn allows for a simple definition of the Chern invariant 



C = 



(19) 



(20) 



It can now easily be seen that Eq. ( 18 1 is invariant with respect to shifts in the energy 



zero (as it must), as long as the Chern invariant vanishes, which we have assumed so 
far. From this, one can heuristically find an expression for the orbital magnetization, 
even for cases in which the Chern invariant is not zero, by simply enforcing the 



constant energy shift invariance again. It follows that a suitable extension of Eq. ( 18 1 
is 



M 



orb 



^ImV 
2hc ^ 



d 3 k 
(2tt) 3 



/nk (<9 k Unk| X (Hk + £ nk - 2n) |<9 k U„k) 



(21) 



where fj, is the chemical potential. The factor of 2 in front of p, is necessary to 
compensate the shift of the energy zero, which affects the Hamiltonian as well as 
the eigenvalues. 



Equation (21) together with (20) lead to 



dpi 



ch{2n) 2 



(22) 



for any insulator and [i in the gap, suggesting that the magnetization changes 
linearly with /i if the Chern invariant is non-zero. This peculiar result is a remnant of 
the surface states present in Chern insulators^ and has been proven numerically.^ 
Without proof, but based on numerical simulations, it can also be shown that 



Eq. (21) holds for metallic systems,^ and is thus a generally valid expression for 



all solids. Its correctness is further supported by alternative derivations, which are 
presented in the following sections. 
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3.4. Semi- classical derivation of the orbital magnetization 

Both of the previous sections derived a fully quantum-mechanical theory of the or- 
bital magnetization in solids based on a transformation to well-localized Wannicr 
functions to circumvent the problems of the position operator in extended systems. 
At the same time the Wannier function approach was developed, a separate research 
group — completely independent — derived an identical expression of the orbital mag- 
netization through very elegant arguments concerning the semi-classical equations 



of motion for Bloch electrons, as described in detail in Ref. 75 . Its derivative, 
presented below, and the derivations from the previous sections nicely complement 
each other, each offering insight into certain aspects of the same physics that the 
other cannot provide. 

The derivation starts from the semi-classical equations of motion in Gaussian 
units for an electron wave packet in band n, i.e. 

. 1 dE nk 



h <9k 



k x n nk , (23) 



Sk= -eE(r) - ~rx B(r) , (24) 
c 

where E(r) and B(r) are the external electric and magnetic field, respectively. For 
simplicity, in the following we will assume spatially constant external fields. It is 



interesting to note that the second term on the right hand side of Eq. ( 23 ) is often 
overlooked in elementary textbooks. This "anomalous velocity" containing the Berry 
curvature is important for properties such as the anomalous Hall conductivity.^ 

Closely linked to these semi-classical equations of motion is the phase-space 
volume element AV = ArAk, which Liouville's theorem predicts to be conserved, 
as is schematically depicted in Fig. [2j If we consider the equation of motion for the 



volume element, 123 (l/AV)(dAV/dt) = V r r + V k k, and insert Eqs. Q and (|24| 
we find for the time evolution of the volume element 

AV = AV " - . (25) 

1 + eB • fl nk /hc v ' 

Since tt n k depends on k, it appears that the volume element changes during the time 
evolution of the state variables (r, k), violating Liouville's theorem. The situation 
can be remedied by introducing a modified density of states 

Dnk = tAv(1 + eB • fl nk /hc) , (26) 

such that the number of states in the volume element D^AV is constant again. 
Then, the expectation value of a physical observable O can be written as 

n J 

For systems with time-reversal and inversion symmetry in the spatial wave function 



the Berry curvature is zero and Eq. (27) does not reveal any new physics. But 



in systems where the Berry curvature is non-zero such as crystals with broken 
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k 




Fig. 2. Schematic depiction of the volume element AV = 
ArAk in phase space (r, k). The volume element might 
change its shape as a function of time, but the volume 
itself is conserved according to Liouville's theorem. (Figure 
modeled after Ref. [78] ). 



time-reversal or inversion symmetry, Eq. (27) provides an elegant way to study 
Berry-phase effects. 



In the original paper, Ref. 75 , the authors give several examples of how Eq. ( 27 1 



can be used to derive, in a simple manner, expressions for the electron density or the 
anomalous Hall conductivity. In particular, it can be used to derive an expression 
for the orbital magnetization. To that end, we describe an electron semi-classically 
through a wave packet in Bloch band n. Such a wave packet is found to rotate about 
its center of massp^ leading to a magnetic moment]^] 



m„ k = - — (<9 k U nk | X (if k - £ n k)|#kUnk) 



(28) 



In a weak magnetic field B, the band structure -E nk = (V'nkl-HlV'nk) changes to 



E n k — m„ k • B. If we then express the total energy according to Eq. (27) as 
d 3 k 



E 



E 



/d k 
fhk (1 + eB • fl nk /hc) (E nk - m nk • B) , 



(29) 



the orbital magnetization follows by the negative derivative of the energy with 
respect to B in the zero field limit as 



e 
2hc 



d 3 k 



/nk (dkUnk\ x (i? k + E nk - 2/i) |<9 k W„k) 



(30) 



in perfect agreement with Eqs. (14) and (21) 



Since this semi-classical derivation does not rely on the localization properties of 
Wannier functions, it is valid for all insulators and metals alike, further supporting 



the heuristic arguments from Section 3.3 



d Note that there is a factor of —1/2 missing in the original equation published in Ref. [79] , as 
pointed out in Ref. [7*5] . 
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3.5. Finite-temperature derivation of the orbital magnetization 
and generalization to interacting systems 

The derivatives of the orbital magnetization formula from the previous sections had 
limitations in that they were either only valid for insulators or they were based 
on semi-classical arguments. Here, we review a completely general derivation of 
the same expression, which is free of these limitations and was first presented in 



Ref. 80 



The starting point for this derivation is the thermodynamic definition of the 
orbital magnetization 

where fig = E — TS — uN is the grand canonical potential, B is the magnetic field, 
and Vq is the unit-cell volume. The derivative is to be taken at constant temperature 
T and chemical potential u. It is convenient to first focus on the related quantity 

with K = E — uN. The original orbital magnetization can then be recovered by 
integration of 

M orb = |^M orb , (33) 

where f3 = 1/kT. 

Next, we apply an external, slowly varying field B(r) with vector potential A(r) 
and wave vector q, pointing in the z-direction according to 

B(r) = Bcos(qy)i , A(r) = -B^-^x. (34) 

At the end of this derivation, we will take the limit q 0. If we consider the 
original system to be described by the Hamiltonian Ho with Bloch eigenfunctions 
ip n k and eigenvalues i?„k, we can treat the magnetic field as a perturbation, e.g. 
H = Hq + Hb , with 

ffu = ^[v-A(r)+A(r).v] . (35) 
We can then define the grand-canonical ensemble energy density as 

K ^ = E J (2^)3 /« k Re K V„k(r)J , (36) 

where K is now the operator K = H — uN and K = H — uN. Perturbation theory 
yields in first order a change in K according to 

v—v f d 3 k r 

&K{v) = J2 7^)3 Re [ A ^«k i>rik K IPnk + fnk Ck H B *Pnk 

+ /„k (AV4k K Q V>„ k + < k K o A^ lk ) ] • (37) 
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However, in linear order in B the change in K can also be written as 



AK(r) 



B(r), 



(38) 



such that, through comparison of Eqs. (37) and (38), on e ca n find an explicit ex 



pression for M orD . Note that the first two terms of Eq. (37) do not contribute to 



M or b Inserting the first order change of the wave function ip n k in Eq. (37) and 



taking the limit q — » 0, one finds after some straightforward but tedious algebra^ 

d 3 k 



M, 



orb 



e 
2hc 



(2tt) 3 



/nk (3 k Mnk| X (H k + E nk - 2fJ,)\dk.Unk) 



df n 



9E„ 



:(Enk - |U)(9kU„k| X {E n k - #k)|<9 k Mnk) 



(39) 



For T = 0, the derivative 



g/ak 



becomes a (5-function of (_E„k — /i) and the 
second term in Eq. (39) vanishes, such that the remainder is in perfect agreement 



with Eqs. (14), (21), and (30) from the previous sections. A finite temperature 



expression can be simply derived through integration of Eq. ( 33 ) and yields 




(40) 



where the Berry curvature fl n k and m„k are defined in Eqs. (19) and (28), respec- 
tively. 

In Ref. [80] the authors go on to show that Eq. (40), besides being valid for 



insulators with or without Chern invariants and metals at zero or finite temper- 
atures, is also valid for weak and strong magnetic fields. It further is applicable 
to interacting systems, provided that the single particle energies and wave func- 
tions are computed using current and spin density functional theory (CSDFT)P^ 



As such, Eq. (40) presents the most general derivation and concludes our review of 



derivations for the theory of orbital magnetization. 



4. Practical Aspects of Calculating the Orbital Magnetization 
4.1. k-space derivatives 



The expressions for the orbital magnetization in Eqs. (14), (21), (30), and (40) 
contain the fc-space gradient of the cell-periodic part of the Bloch function, i.e. 
|<9 k w n k) with 3k = 9/9k, the evaluation of which is not necessarily trivial. In the 
following, we review three practical options to compute this derivative. 

The first option can be considered analytical and is based on standard pertur- 
bation theory. The change of the states in first order to a k • p perturbation can be 

e Note that the cross product X is missing in the original expression in Eq. (13) of Ref. [80] . 



Theory of Orbital Magnetization in Solids 15 



expressed as a sum over states, i.e. 

|5 k W„ k ) = 2^ \Umk)-^ — —b — , (41) 

m^tn 

where the velocity operator v is defined in Eq. . The disadvantage of this expres- 
sion is that all states have to be included, even unoccupied ones. However, often 
first-principles computer codes — to save compute time — do not calculate many un- 
occupied states. Also, when plane waves are used, the sum converges very slowly and 



becomes unpractical. Equation (41) is, however, useful in the context of localized 



basis sets or tight-binding models, where the number of states is small. In particular, 



this method has been used to calculate the results presented in Section 4.3 

The second option is numerical in nature, making use of the covariant deriva- 
tive J82|lM] rp^g numer i ca i derivative is not trivial, as can be seen from the following 
example: let q be a vector that connects k with a nearby point in direction a. In 
principle, we would then like to approximate the desired derivative as the finite 
difference 

|<9 Q U„ k ) ~ (l U «k+q) - Ktk-q)) , (42) 

and the expression becomes exact in the limit of |q| — > 0. However, numerically 
this is not a well-defined quantity. While the states |u n k+ q ) and |w„k- q ) can easily 
be obtained through diagonalization of the Hamiltonian at these points, they in 
general will have arbitrary and unrelated overall phase factors, depending on the 
diagonalization routine used. Thus, simply numerically subtracting the states yields 
unpredictable results. 

The problem of the phase can be circumvented by replacing the normal derivative 
with the covariant derivative^!!!] 

\d a u nk ) -> \d a u n u) = Qk|<9 a w„ k ) • (43) 

Here, Qk projects onto the empty states at point k and is defined as the com- 
plement of Vic, i.e. Qk = 1 — T^k, which projects onto the occupied manifold 
7\ = Y^n° C l w nk)( u nk|- Note that the definition of a derivative on such a mani- 
fold is not unique and the covariant derivative is a choice; other choices exist, which 
also have components in the occupied states. We can then define dual states |Mnk+q) 
as a linear combination of the occupied states |M m k+q) at point k + q as 

|w»k+q) = X^k.k+qWK'k+q) , (44) 
rn 

where S is the overlap matrix 

(<Sk,k+q) > • (45) 

For sufficiently small q, S becomes nearly diagonal, containing only the arbitrary 
overall phase factors of (u n k| and |w n k+q) in its diagonal elements. The covariant 
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derivative can then simply be written as 

\d a U nk ) PS ^r— (j£t„k+q) - |£i„k-q)) , (46) 

which is manifestly gauge independent in the sense that it does not depend on 



the choice of phases for the |u„k± q )- It thus resolves the phase issue of Eq. (42 1 
and is numerically accessible. Note that the covariant derivative \d a u nk ) still has 
the same overall phase factor as |u„k), which then cancels when evaluating matrix 
elements such as the orbital magnetization. The effect that the overlap matrix plays 
in canceling the arbitrary overall phase can be seen more clearly if only one band 



is considered, as is done explicitly in Eq. (26) of Ref. 82 



The dual states enjoy the property of being orthogonal to the original states 

(w„k|Umk±q) = S nm , (47) 

from which follows immediately that 

(UnklBaUmk) = — r (S nm - 8 nm ) = . (48) 

^19.1 

In other words, the covariant derivative is orthogonal to the occupied states, as 



evident by the projector Qk in Eq. (43) 



The covariant derivative requires that the number of occupied states at k and 
k ± q is the same, i.e. the material has to be an insulator. This is a remnant of 
the fact that it was first developed in the context of the electric polarizationjUHMl 
which is only defined in insulators. However, it is conceivable that the formalism of 
the covariant derivative can be extended to metallic systems. 

The third option for calculating fc-space derivatives is based on linear-response 



theory. If we multiply Eq. (41 ) from the left with — E n k, we find 

(flk - E nk ) |<9kw„k) = -Q«kv|n„ k ) , (49) 
where Q„k = 1 — \u nk }(u n k\- This is a Sternheimer equation, similar to Eq. (25) 



in Ref. 84 . The advantage of this expression is that the expensive sum over all 
states vanishes and the right hand side can be evaluated with only knowledge of 
the occupied states. The equation then forms a linear system, which can easily 
be inverted to find \d k u nk ) by iterative algorithms, such as a conjugate-gradient 
approach. This linear-response approach can be used to calculate the orbital mag- 
netization in real, metallic systems and, in particular, has been applied to calculate 
the orbital magnetization in Fe, Co, and Ni, discussed in Section [5] 



4.2. Calculating the orbital magnetization 
in a pseudopotential context 

Many first-principles electronic-structure codes use plane waves as basis functions 
for their many conveniences. However, to allow for a reasonable basis-set size the 
true Coulombic potential of the ions is usually replaced by an effective potential, i.e. 
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pseudopotential. Since the expressions for the orbital magnetization in Eqs. ( 14 1 



(21), (30), and (40) are strictly speaking only valid in an all-electron picture, they 



have to be modified to be used within a pseudopotential framework. The corre- 
sponding theory has been recently developed for the special case of norm-conserving 
pseudopotentials in Ref. [8], and we present here a summary thereof. 

In the pseudopotential formalism, the pseudo Hamiltonians H and act on 
the pseudo wave functions \ip n k) and |M n k) with eigenvalues E n ^. As a first approx- 



imation, either of the expressions in Eqs. (14), (21), (30), or (40) can simply be 
evaluated using the pseudo quantities H^, |tt n k), and E^ instead of the all-electron 
ones H, |u„k), and E n ^. This results — as we shall see in Section [5] — in only a small 
error in the orbital magnetization, on the order of a few percent. 

The true wave functions can be recovered from the pseudo ones by means of 
the projector augmented waves (PAW) transformation |V>nk) = T|V>nk)^ I n the 
presence of magnetic fields, the PAW transformation has to be modified, resulting 
in a gauge including projector augmented waves (GIPAW^P formalism with the 
corresponding transformation |V>nk) = TB|V>nk)f^ The pseudo Hamiltonian then is 
given through H = T%BT B . If we define the orbital magnetization as the derivative 
of the total energyj^] it follows that 

1 dE tot 



Vn dB 



B=0 



E 



i 

-E 



1 

Vn 



E 



d 3 k 
(2^r)3 

d 3 k 
(2tt) 3 

d 3 k 
(2^)3 



B=0 



d 

/nk ^(V'nkl.HlV'nk) 

d - 

/nk ^(V , nk|-H r |V'nk) B=0 



/nk (i>nk\^\i>nk) 



B=0 



(50) 



where we have used the Hellmann-Feynman theorem in the last step. Explicitly 
calculating the pseudo Hamiltonian through the GIPAW transformation, taking 
the B derivative, and finally taking the B = limit yields 



wi 



itrH 



AM bare 
orb 



AM para 



M orb = M^ c + AM^ C + AM^ a + AM^b , 



E / 7?T3 /nk <^k|(R - r) x i [r - R, V^} |u nk ) , (52) 



(51) 



2hcV Q 



9'e 
16m^c 3 Vo 



E 

nR 



AM^ = - 



g'He 



E 



16m 2 e c 3 V ^ J (2tt) 3 

n R, 



d 3 k 
(2^)3 

d 3 k 



/nk (u„ k |(R - r) x i [r - R,F£ L ] \u nk ) , (53) 

/nk^nklE^l^k) • (54) 



Here, we are not including the trivial spin-Zeeman term in the Hamiltonian; if we did, the energy 
derivative would correspond to the total magnetization. 
g Note that there is a factor of — 1/Vb missing in Eqs. (10)-(13) in Ref. [jj]. 
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Here, g' « 2.004639,^ m e is the electron mass, is the non-local part of the 
pseudopotential in separable form, and and E^ L are the separable non-local 
GIPAW projectors at atomic site R. If we denote the GIPAW projector as \p n R) 
and the all-electron and pseudo partial waves as |^> n R.) an d \<p n R), respectively, the 
latter terms can be written as 



nm 



(vW|W x p\c/) mR ) - (<f> n n\VV x p\<p mR ) (p mR \ , (55) 
a|r x (cr x VV)\(j) mR ) - (<?W|r x (<r x W)|0 mR ) (p mR \ , 



where <x are the Pauli matrices and V and V are the screened atomic all-electron 
potential and the local pseudopotential, respectively. 

The M^ e term on the right hand side of Eq. (51 ) is indeed the earlier derived 
expression for the orbital magnetization — Eqs. (14), (21 1, (30 1 — simply evaluated 
for the pseudo Hamiltonian and pseudo wave functions |it ra k) with eigenvalues 
E n u- This "bare" term is the largest contribution to the orbital magnetization and 
the other terms in Eqs. ( 52 1— ( 54 ) can be considered corrections to the bare value. 
At first sight it might appear that these correction terms are plagued from the 
same ill-defincdness of the position operator r in periodic systems, as described in 
Section 3.1 But, the calculation of these corrections is actually trivial — even in 
extended systems — since the action of the non-local operators , , and E^ L 
is non-zero only inside spherical regions around the atom at R. 

First calculations of the orbital magnetization for Fe, Co, and Ni have been 
carried out in Ref. |8; using this pseudopotential approach. We report results for 
these calculations in Section [5j where Table [l] shows explicit results for all terms 
that contribute to the orbital magnetization in Eq. (51). 



4.3. Single k-point derivation of the orbital magnetization 



The expression for the orbital magnetization in Eq. (21) is a fc-space integral over 
the Brillouin zone. The size of the Brillouin zone is inversely related to the size of 
the unit cell in real space. As such, it is interesting to study how large a supercell 
has to be such that the Brillouin-zone integral is reasonably well approximated by 
a single fc-point. Assuming a large supercell, in Ref. 89 an effective single fc-point 
expression is derived, which is particularly useful since most simulations for non- 
crystalline systems, including Car-Parrinello simulations,^ are routinely performed 
using a large supercell in combination with only one fc-point. We present a simpli- 
fied derivation here and refer the reader to the original manuscript in Ref. 89 for 
further details. 



The starting point is Eq. (21), slightly rewritten to include the antisymmetric 
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tensor e 7Q/ g and the shortcut d a = d/dk a ^ 

v—v e f d 3 k 

M orb , 7 = 2J e 7«/3 ^ Im / T^pjl /™ k (<9aU«k|^k + E nk - 2ii\dpu nk ) . (56) 

For a sufficiently large supercell of volume V, this expression can be approximated 
by replacing the fc-space integral by the single fc-point k = 0, 

M olhn k, ^ € ja/3 -^-^ Im f n0 (d a u n0 \H + E n0 ~ 2n\d f3 u n0 ) . (57) 

This approximation is possible because the integrand for the orbital magnetization 
itself is gauge invariant.^ Note that the analog case of the electric polarization, 
where only the integral as a whole is gauge invariant, is more complicated.^ If we 
denote the shortest reciprocal lattice vectors of the supercell as hj and the derivative 



in the same directions as dj, Eq. (57) can be written as 

M orb « ^2 tiji bi|bj-||b|| 6 - 3 Im f n0 (djU n0 \H + E n0 - 2fi\diu nQ ) , (58) 
ijl,n ^ ' 

where the fc-space derivatives of |m„o) can be understood as the finite difference 

\diu n0 ) = lim -T-^-r(|u n Abi) - \u n o)) ■ (59) 

Such derivatives are numerically difficult due to gauge issues and special care has 
to be taken as described in Section [4.1| Approximating this expression with A = 1 



and using the dual states defined in Eq. (44), the orbital magnetization becomes 
M orb ~ ^ €iji b t 6 Im /„ (unhj \H + E n0 - 2(i\u nbl ) . (60) 

ijl,n ' 

Note that the matrix elements (u„o| ■ • • \u n o) and (u n b - \ ■ • ■ \uno), due to the orthog- 



onality of the dual states in Eq. ( 47 1 , are purely real and vanish after taking the 
imaginary part. 

Since the periodicity in /c-space is described by the vectors bj, the next nearest 
neighbor points necessary for these derivatives can all be related to the k = point 
via jwnbj) = e~ 3 " r |it n o). The resulting expression for the orbital magnetization is 
thus indeed a single fc-point formula in the sense that the Hamiltonian has to be 
diagonalized only once at k = to obtain the eigenvalues E n Q and the states |u n o)j 
the orbital magnetization can then be computed only based on these quantities. 

The efficiency of this single fc-point formula has been tested in the context of 
a simple two-dimensional tight-binding model developed by HaldaneP^ This model 
has no macroscopic magnetization and breaks time-reversal symmetry through a 
complex second-nearest neighbor hopping, resulting in staggered magnetic fluxes 
that cancel throughout the unit cell. The convergence of the orbital magnetization 
with respect to the size of the supercell can be seen in Fig. [3] 



h The corresponding expression, Eq. (2) in Ref. [89], has an erroneous minus sign, which propagates 
through the remainder of that reference. Also, several sums over n are missing in Ref. [89] . 



20 T. Thonhauser 



0.082 




0.074 1 1 1 1 1 

0.05 0.1 0.15 0.2 

1/L 

Fig. 3. Orbital magnetization of the Haldanc 
units of the Bohr magneton fiB as a function of the su- 
percell size L. The exact value is compared with the re- 
sults obtained by evaluating the single fc-point formula from 
Eq. \<o0\ , using the numerical and analytical methods for 
the fc-space derivatives outlined in Section |4.1| (Reprinted 
with permission from Ref. 89 ; © 2007 American Physical 
Society). 



4.4. Wannier interpolation of the orbital magnetization 

An efficient, alternative numerical approach to evaluate the orbital magnetization 
that targets metallic materials specifically is currently being developed.^ The basic 
idea is to first perform an ab initio ground-state calculation on a relatively coarse 
fc-mesh that yields ground-state wave functions and energies. From these, partially 
occupied disentangled Wannier functions are constructed,^ which then serve as a 
compact tight-binding-like parameterization of the ab initio results. The Wannier 
representation can then be used to interpolate quantities of interest on (almost) 
arbitrarily fine /c-meshes and evaluate the corresponding Brillouin-zone integrals 
easily. This method of Wannier interpolation has been developed in the context of 
the anomalous Hall effect in ferromagnets, where A:-meshes of up to 300 x 300 x 300 
grid points were necessary to converge the Brillouin-zone integral for the anoma- 
lous Hall conductivity^ — a quantity which has much in common with the orbital 
magnetization. Although the Brillouin-zone integrand of the orbital magnetization 
in Eqs. (14), (21 1, (30), and (40) is not as spiky as the one for the anomalous Hall 
conductivity, a Wannier interpolation scheme still provides advantages. At the time 
of writing, the corresponding manuscript discussing the interpolation of the orbital 
magnetization in Ref. [92| is under review, and we refer the reader to the published 
version — once available — for discussion, results, and further details. 
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4.5. DFT specific aspects of calculating the orbital magnetization 

The formalisms collected in this review article describe the calculation of the orbital 
magnetization in periodic solids. As such, they can be used in conjunction with 
any first-principles method suitable for solids. While other successful approaches 
exist, density functional theory (DFT^S has been the primary tool for modeling 
many molecules and solids}^ in fact, the vast majority of all electronic-structure 
calculations today are DFT It is thus appropriate to consider some aspects of the 
calculation of the orbital magnetization that are DFT specific. 

One of the most important "ingredients" for DFT calculations is the exchange- 
correlation functional used. As such, it is interesting to study which exchange- 
correlation functional gives results for the orbital magnetization in best agreement 
with experiment. However, since the theory of orbital magnetization in solids is 
rather young, a systematic study of the orbital magnetization's dependence on the 
functional has not been performed yet. The results for Fe, Co, and Ni discussed in 
the next section suggest a slightly better performance of PBEp3 compared to LDAj^l 
but that might also be related to PBE's better overall performance in solids. 

However, even if we had access to the true exchange-correlation functional, DFT 
is only guaranteed to produce the correct ground-state density; a similar claim for 
the orbital currents cannot be made. Thus, in principle, the calculation of the orbital 
magnetization — which is defined in terms of orbital currents — within DFT is not 
guaranteed to yield correct results. It seems that a formulation that also includes 
currents, such as current density functional theory (CDFT)J^3 should systematically 
improve results for the orbital magnetization. But, it remains to be seen if appropri- 
ate current functionals, in conjunction with CDFT, do indeed yield improved results 
over standard DFT. Alternatively, time-dependent current density functional the- 
ory (TD-CDFT^HEnn] or the better developed time-dependent density functional 
theory (TD-DFTjP^ could in principle be employed. But again, at this point it 
is not clear that they provide any practical advantage for calculating the orbital 
magnetization over standard DFT. A somewhat more detailed discussion of these 



aspects is given in Ref. 73 



5. Orbital Magnetization in Fe, Co, and Ni 

Calculations of the orbital magnetization using the full theory have first been re- 
ported in Ref. [8]. In that paper, the authors calculate the orbital magnetization of 



Fe, Co, and Ni within DFT by evaluating Eq. (21 1 in the pseudopotential context of 



Section 4.2 using the linear-response method for the fc-space derivative described 



in Section 4.1 This framework was implemented in the first-principles electronic- 
structure code PWSCF, which is part of the Quantum-Espresso packageP^H At 
the time of writing, this implementation is not yet publicly available. 

Results for the orbital magnetization and its split-up into the bare and correc- 
tion terms according to Eqs. (51)-(54) are collected in Table [lj values for Af or b 



are taken from Ref. 18 and values for the bare and correction terms are taken 
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Table 1. Orbital magnetization of Fe, Co, and Ni in units of the Bohr magneton fig per atom 
parallel to the spin, for different spin directions and different exchange-correlation (XC) functionals. 
The total orbital magnetization is given by M orb = M^ c + AM^ C + AM' r " a + AM^, 
where the contributions are defined in Eqs. ( ] 5 1 [ > — ( |54p . Values are taken from Ref. [8] and its 
supplemental materials; sec footnote in main text. Experimental values are taken from Rcfs. 3 
and [59| . (Reprinted with permission from Ref. [8]; © 2010 American Physical Society). 



Metal 


dir. 


exp. 


XC 


"^orb 


M barc 
orb 


AM bare 
orb 


AM p r 

orb 


AM dia 
orb 




[001]* 


0.081 


PBE 


0.0658 


0.0639 


0.0000 


0.0016 


0.0003 


fecc-Fe 


LDA 


0.0642 


0.0616 


0.0005 


0.0017 


0.0004 


[111] 




PBE 
LDA 


0.0660 
0.0633 


0.0637 
0.0609 


0.0005 
0.0005 


0.0015 
0.0015 


0.0003 
0.0004 




[111]* 


0.120 


PBE 


0.0756 


0.0722 


0.0018 


0.0014 


0.0002 


fcc-Co 


LDA 


0.0741 


0.0706 


0.0019 


0.0014 


0.0002 


[001] 




PBE 
LDA 


0.0660 
0.0642 


0.0629 
0.0608 


0.0016 
0.0016 


0.0013 
0.0016 


0.0002 
0.0002 




[001]* 


0.133 


PBE 


0.0957 


0.0908 


0.0032 


0.0014 


0.0003 


hep-Co 


LDA 


0.0924 


0.0875 


0.0032 


0.0014 


0.0003 


[100] 




PBE 
LDA 


0.0867 
0.0837 


0.0822 
0.0792 


0.0028 
0.0029 


0.0015 
0.0013 


0.0003 
0.0003 




[111]* 


0.053 


PBE 


0.0519 


0.0494 


0.0017 


0.0007 


0.0001 


fec-m 


LDA 


0.0545 


0.0519 


0.0019 


0.0007 


0.0000 


[001] 




PBE 
LDA 


0.0556 
0.0533 


0.0527 
0.0505 


0.0022 
0.0020 


0.0006 
0.0007 


0.0001 
0.0001 



* denotes the experimental easy axis. 



from the supplemental materials of the same reference^ The results clearly show 
that, as pointed out in Section [2j the orbital magnetization in Fe, Co, and Ni is 
indeed a small effect of only a few percent compared to their spin magnetization 
of 2.083 hb (Fe), 1.523 [i B (Co), and 0.518 [i B (Ni)P^It is also interesting to see 
that the orbital magnetization is well described by the bare term alone and 

the pseudopotential corrections AM^™, AM^ a b ra , and AM*J are usually small. 
Although only a small effect, it can also be seen from these results that PBE gives 
slightly better results than LDA. However, that might also be related to the better 
performance of PBE in solids in general. 

With the full theory of orbital magnetization available, we can now assess the ac- 
curacy of the often-used muffin-tin approximation, mentioned in Section [3. 1| In this 
approximation, the orbital contribution can be simply calculated using Eq. @, in- 
tegrating over non-overlapping spheres centered around the atoms. Since the sphere 
volumes are finite systems, the problem of the ill-defined position operator is cir- 



1 Reference [5] and its supplemental materials only contain PBE results. The values for LDA results 
were obtained from the authors of Ref. |8| through personal communication. 
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Table 2. Orbital magnetization of Fe, Co, and Ni in fig /atom parallel 
to the spin, for different spin directions. The total orbital magnetization 
is here split up into contributions coming from the muffin-tin spheres, i.e. 
Eq. \61) , and the interstitial. Calculated values are taken from Ref. [8| and 
experimental values are from Refs. [3] and [59| . (Reprinted with permission 
from Ref. 18]; © 20f0 American Physical Society). 



Metal dir. exp. ^oib interstitial muffin tin 



bcc-Fe 
fcc-Co 
hep-Co 
/ec-Ni 



[001]* 

[111] 

[111]* 

[001] 

[001]* 

[100] 

[111]* 

[001] 



0.081 
0.120 
0.133 
0.053 



0.0658 
0.0660 
0.0756 
0.0660 
0.0957 
0.0867 
0.0519 
0.0556 



0.0225 
0.0216 
0.0122 
0.0064 
0.0089 
0.0068 
0.0008 
0.0047 



0.0433 
0.0444 
0.0634 
0.0596 
0.0868 
0.0799 
0.0511 
0.0509 



* denotes the experimental easy axis. 



cumvented. In particular, one can calculate 
« = ~ E / (03 U M* x v| Unk ) MT (61) 



7^3 / d\ /„ k < k (r) [r x (-ihV + Kk)/m e ] u nk (r) , 
\ Z7T ) Jmt 

where we integrate inside all muffin-tin (MT) spheres in real space. This muffin-tin 
contribution was calculated in Ref. [8] for a sphere radius of 2 Bohr. The interstitial 
contribution is then defined as the total orbital magnetization minus the muffin-tin 
contribution. Results for Fe, Co, and Ni are collected in Table [2j Surprisingly, the 
results show that up to 34% of the orbital magnetization in Fe, up to 16% in Co, and 
still up to 8% in Ni originate from the interstitial. As such, it is not surprising that 
calculations of the orbital magnetization — using the muffin-tin approximation — 
strongly underestimate the effect HHHS These results, a posteriori, provide indeed a 
strong justification for developing the complete theory of orbital magnetization. 

The overall results for the orbital magnetization in Table [2] have also been con- 
firmed with an independent implementation of the Wannier interpolation method 
presented in Section [4~4| and Ref. 92 . On the other hand, preliminary calculations 
by another group, using an implementation in the Wien2k codep^ disagree with 
the results in Table [2J however, those results were never published?^ 



6. Applications of the Orbital Magnetization 

While the orbital magnetization of materials itself is a useful concept, its change 
is of much practical interest, as it can be related to several experimental probes. 
The calculation of changes of the orbital magnetization has been possible for about 
15 years now, using linear-response methods. 52 56 But, now that the orbital mag- 
netization itself can be calculated, numerous possibilities open up to develop new 
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approaches — much simpler than linear-response methods — using finite differences. 
Such approaches might prove particularly useful in situations where linear-response 
calculations are cumbersome or impossible, or in conjunction with more complex 
methods of treating exchange and correlation effects, such as hybrid functionals, 
DFT+U, exact exchange, or beyond-DFT approaches. In the following sections we 
focus on calculating the nuclear magnetic resonance shielding tensor and the electron 
paramagnetic resonance g-tensor as a finite difference of the orbital magnetization. 



6.1. Nuclear magnetic resonance as a derivative 
of the orbital magnetization 

Nuclear magnetic resonance (NMR) is one of the most important experimental tech- 
niques used to determine the structure of molecules, liquids, and other disordered 
systems; it has thus evolved into one of the most widely used methods in structural 
chemistryllHHHD Similar to the orbital magnetization itself, methods to calculate the 
NMR response for molecules and clusters, i.e. finite systems, were developed early 
on in the quantum-chemistry communityp^ but calculations for periodic solids were 
impossible, as the inclusion of a constant magnetic field requires a vector potential 
that breaks translational symmetry. While linear-response frameworks to calculate 
the NMR shielding in solids have been available for a while JH1HMH] knowledge of the 
orbital magnetization itself provides a much simpler way to calculate the shielding, 
referred to as the converse NMR approach. The corresponding formalism is de- 
scribed in detail in Ref. [7]. 

When a sample is put in an external magnetic field B cxt , a field B lnd is induced, 
resulting in a total magnetic field of B* ot = Bf d + B cxt at the site of atom i. 
The NMR shielding tensor <Tj for atom i is then defined as B" ld = — cr s: • B oxt , or 
equivalently as 

dB ind 

= -Q^k ■ ( 62 ) 

We now consider an artificial magnetic dipole at site i. Such a dipole would have 
the energy E = m i ■ B* ot , which leads to 

B*°\ = -dE/dm aii . (63) 



Starting from Eq. (62), we can now use a thermodynamic relationship between 



mixed partial derivatives and insert Eq. ( 63 ) to find 



~ B ° a ) _ x d dE 

a «P,i - <W* 



<95° xt ap dBf*dm a>i 

= + dm a>i BBf = 6afi ~ V °^- ' (64) 

where Vq is the unit-cell volume and we have written the macroscopic orbital mag- 
netization as M or b„g = —(l/Vo)dE/dB t p t . By doing so, in the last step we have 
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made the assumption that the external held can be replaced by the total held, which 
corresponds to a particular choice of geometry for the sample shape. However, the 
results for other sample shapes can easily be recovered through knowledge of the 
susceptibility.^ 



Equation (64) expresses the shielding tensor in terms of changes in the orbital 
magnetization due to the presence of the artihcial dipole m.; , and hnite differences 
can readily be calculated. To that end, the dipole needs to be included in the 
system's Hamiltonian via its vector potential, for which a suitable form has been 
derived^ — a minimal change to existing codes, compared to a linear-response frame- 
work. Note that this approach circumvents the difficulties related to including hnite 



external magnetic helds and the choice of gauge origin. In order to evaluate Eq. ( 64 I 



in a pseudopotential framework, a GIPAW transformation similar to Section 4.2 has 



to be performed, which includes the vector potential of the artihcial dipole in the 



Hamiltonian; the corresponding framework has been derived in Ref. 109 . 

The converse NMR approach has been implemented in the first-principles 
electronic-structure code PWSCF, which is part of the Quantum-Espresso pack- 
a g e |i02j j^. ^ e tj me of writing, this implementation is not yet publicly available. 



Calculations have been performed within DFT by evaluating Eq. (64) through fi- 
nite differences in the GIPAW pseudopotential context using the linear-response 
method for the fc-space derivative described in Section |4.1| This implementation 
has already been successfully applied to calculate H, C, F, O, P, and Si NMR shifts 
in a variety of systems such as small molecules, polycyclic aromatic hydrocarbons, 
bulk water, and selected Si0 2 crystalsIOElEin] 

6.2. Electron paramagnetic resonance as a 
change of the orbital magnetization 

Electron paramagnetic resonance (EPR) is in nature similar to NMR, except that 
the focus is on electronic spins instead of nucleic spins. Most stable molecules have 
no unpaired electrons, such that EPR is limited to paramagnetic materials and is 
thus less widely used. The property of interest is then the electron g-tensor. As in 
the case of NMR, a linear-response framework exists to calculate the g-tensorP^l 
However, knowledge of the orbital magnetization allows for a direct calculation.^ 

The g-tensor deviation Ag a p from the free-electron value of g e — 2.002 319 can 
be calculated through a change in the orbital magnetization aijSKSl 

1 d 

Ag Q/ 3 = ^-M or b, , (65) 

where S is the total spin. This approach has been used to study the microscopic 
structure of radicals and paramagnetic defects in solidsPln particular, calculations 



have been performed at the DFT level by evaluating Eq. ( 65 1 through hnite differ- 
ences according to a spin hip in the GIPAW pseudopotential context described in 
Section [4.2| using the linear-response method for the fc-space derivative described in 
Section |4~T| It is interesting to note that, while for the orbital magnetization the 
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pseudopotential correction terms in Eqs. (52)-(54) play only a minor role (see Ta- 
ble [l] on page 22), the correction terms play a more significant role when calculating 
A# Q/3 . 



6.3. Other derivatives of the orbital magnetization 



Derivatives of the orbital magnetization also play an important role in other areas 
of condensed matter physics. We conclude this section with a non-exhaustive list, 
mentioning several recent examples. 

The magnetic susceptibility of a material is defined as the derivative of the total 
magnetization with respect to an external magnetic field. As such, the derivative of 
the orbital magnetization contributes to the magnetic susceptibility in general.^ 



Furthermore, in Ref. 12 the author links the spin-Hall conductivity in insulators 



through a Stfeda-like formula to the magnetic susceptibility, and in particular, to 
the derivative of the orbital magnetization with respect to a magnetic field. Finite 
differences calculations of the orbital magnetization were carried out for a simple 



two-band model of graphene evaluating Eq. (21 ). 

As pointed out in Ref. [13 , it might be possible to identify non-abelian quantum 
Hall states by experimental measurements of the temperature dependence of the 
orbital magnetization. This is of particular interest in view of the recently developed 



analytical formula for the temperature dependent orbital magnetization in Eq. ( 40 ) 



In insulators with broken time-reversal and inversion symmetry, in first order an 
electric field E can induce a magnetization M and a magnetic field B can induce 
an electric polarization P. The coupling is described by the linear magnetoelectric 
polarizability a a p p3 defined as 



B=0 



dB p 



(66) 



E=0 



Orbital contributions to these derivatives are discussed in Refs. [9] and 10 



7. Conclusions and Outlook 

The theory of orbital magnetization in solids is still in its infancy, being developed 
only in 2005. Its "older brother" — the modern theory of electric polarization — has 
been around for over a decade longer EMUI The modern theory of electric polarization 
was hugely successful, as the number of citations of its fundamental papers shows. 
The theory was a breakthrough because it allows calculations of the polarization 
as a bulk property, without the need for expensive calculations on slabs or clusters. 
It is difficult to predict whether the theory of orbital magnetization will enjoy the 
same level of success and widespread use; in light of the fact that magnetism in 
many materials is well described with the spin contribution only, for which theories 
have existed for a long time, it seems unlikely that it will. 

However, the true success of the modern theory of electric polarization came 
about when its formalism was widely available in first-principles electronic-structure 
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codes. In fact, today its formalism is present in almost all of the popular and widely 
used codesP^ Currently, the situation is different for the theory of orbital magneti- 
zation. At the time of writing, prototype implementations only exist for the PWSCF 
part of the Quantum-Espresso package^ 221 and Wien2kP^1 which are not even 
publicly available yet. Implementations in Vasp AbinitP^ and ADpE^are in a 
development stage. Once more implementations become widely available, the theory 
of orbital magnetization can be expected to play an important role in condensed 
matter physics and materials science. 

Acknowledgements 

This review article was written while on sabbatical at the Department of Materials, 
University of Oxford. I would like to thank Wake Forest University for providing this 
opportunity and Prof. Nicola Marzari for being a wonderful host during that time. 
I am truly indebted to Dr. Davide Ceresoli for countless discussions concerning 
the orbital magnetization, for proofreading this manuscript, and for his everlasting 
patience with my questions. I am further grateful for enlightening discussions with 
Profs. David Vanderbilt, Raffaele Resta, and Qian Niu. Finally, I am also grateful 
to Loah Stevens, who helped with the literature search and organization of the 
references for this review article. 

References 

1. Daniel C. Mattis, The Theory of Magnetism Made Simple: An Introduction to Phys- 
ical Concepts and to Some Useful Mathematical Methods, p. 1 (World Scientific Pub- 
lishing, Singapore, 2006). The first chapter of this book gives a very nice introduction 
to the history of magnetism and is available free of charge through the Wold Scientific 
website at http : / /www . worldscibooks . com/etextbook/5372/5372_chap01 . pdf 

2. Lucretius Carus, De Rerum Natura, 1st century B.C. References are to vv. 906 ff., 
in the translation by Th. Creech, London, 1714. Pliny, quoted in W. Gilbert, De 
Magnete, trans. (Gilbert Club, London, 1900), rev. ed., p. 8 (Basic Books, New 
York, 1958). 

3. A. J. P. Meyer and G. Asch, J. Appl. Phys. 32, S330 (1961). 

4. S. Qiao, A. Kimura, H. Adachi, K. Iori, K. Miyamoto, T. Xie, H. Namatame, M. 
Taniguchi, A. Tanaka, T. Muro, S. Imada, and S. Suga, Phys. Rev. B 70, 134418 
(2004). 

5. H.J. Gotsis and I.I. Mazin, Phys. Rev. B 68, 224427 (2003). 

6. J. W. Taylor, J. A. Duffy, A.M. Bebb, M. R. Lees, L. Bouchenoire, S. D. Brown, and 
M.J. Cooper, Phys. Rev. B 66, 161319(R) (2002). 

7. T. Thonhauser, D. Ceresoli, A. A. Mostofi, N. Marzari, R. Resta, and D. Vanderbilt, 
J. Chem. Phys. 131, 101101 (2009). 

8. D. Ceresoli, U. Gerstmann, A. P. Seitsonen, and F. Mauri, Phys. Rev. B 81, 
060409(R) (2010). 

9. A. Malashevich, I. Souza, S. Coh, and D. Vanderbilt, New J. Phys. 12, 053032 (2010). 

10. A. M. Essin, A. M. Turner, J. E. Moore, and D. Vanderbilt, Phys. Rev. B 81, 205104 
(2010). 

11. A.M. Essin, J.E. Moore, and D. Vanderbilt, Phys. Rev. Lett. 102, 146805 (2009). 



28 T. Thonhauser 



12. S. Murakami, Phys. Rev. Lett. 97, 236805 (2006). 

13. N.R. Cooper and A. Stern, Phys. Rev. Lett. 102, 176807 (2009). 

14. R.D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, R1651 (1993). 

15. D. Vanderbilt and R.D. King-Smith, Phys. Rev. B 48, 4442 (1993). 

16. M. V. Berry, Proc. R. Soc. London, Ser. A 392, 45 (1984). 

17. D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 1959 (2010). 

18. R. Resta, J. Phys.: Condens. Matt. 22, 123201 (2010). 

19. P. Blohski, A. Lehnert, S. Dennler, S. Rusponi, M. Etzkorn, G. Moulas, P. Bencok, 
P. Gambardella, H. Brune, and J. Hafner, Phys. Rev. B 81, 104426 (2010). 

20. B. Lazarovits, L. Szunyogh, and P. Weinberger, J. Magn. Magn. Matt. 240, 331 
(2002). 

21. A. Aldea, V. Moldoveanu, M. Nita, A. Manolescu, V. Gudmundsson, and B. Tanatar, 
Phys. Rev. B 67, 035324 (2003). 

22. H. A. Diirr, S. S. Dhesi, E. Dudzik, D. Knabben, G. van der Laan, J. B. Goedkoop, 
and F.U. Hillebrecht, Phys. Rev. B 59, R701 (1999). 

23. J. Bartolome, M. Garcia, F. Bartolome, F. Luis, R. Lopez-Ruiz, F. Petroff, C. Deran- 
lot, F. Wilhelm, A. Rogalev, P. Bencok, N. B. Brookes, L. Ruiz, and J.M. Gonzalez- 
Calbet, Phys. Rev. B 77, 184420 (2008). 

24. P. Gambardella, J. Phys.: Condens. Matt. 15, S2533 (2003). 

25. J. Hong and R. Q. Wu, Phys. Rev. B 70, 060406(R) (2004). 

26. B. Lazarovits, L. Szunyogh, and P. Weinberger, J. Magn. Magn. Matt. 272-276, 
1658 (2004). 

27. P. Gambardella, A. Dallmeyer, K. Maiti, M. C. Malagoli, W. Eberhardt, K. Kern, 
and C. Carbone, Nature 416, 301 (2002). 

28. T. Oka and H. Aoki, Phys. Rev. B 79, 081406(R) (2009). 

29. J. Liu, Z. Ma, A. R. Wright, and C. Zhang, J. Appl. Phys. 103, 103711 (2008). 

30. D. Xiao, W. Yao, and Q. Niu, Phys. Rev. Lett. 99, 236809 (2007). 

31. W. Yao, D. Xiao, and Q. Niu, Phys. Rev. B 77, 235406 (2008). 

32. D.R. Faulhaber and H. W. Jiang, Phys. Rev. B 72, 233308 (2005). 

33. F. von Oppen, D. Ullmo, and H. U. Baranger, Phys. Rev. B 62, 1935 (2000). 

34. I. Meinel, D. Grundler, T. Hengstmann, C. Heyn, D. Heitmann, W. Wegscheider, 
and M. Bichler, Physica E 6, 731 (2000). 

35. Z. Wang and P. Zhang, Phys. Rev. B 76, 064406 (2007). 

36. Z. Wang, P. Zhang, and J. Shi, Phys. Rev. B 76, 094406 (2007). 

37. F. Cheng, W. Zhi-Gang, L. Shu-Shen, and Z. Ping, Chin. Phys. B 18, 5431 (2009). 

38. M. Todorova, L. M. Sandratskii, and J. Kiibler, Phys. Rev. B 63, 052408 (2001). 

39. M. S. S. Brooks, O. Eriksson, L. Severin, and B. Johansson, Physica B 192, 39 (1993). 

40. J.F. Annett, B. L. Gyorffy, and K.I. Wysokihski, New J. Phys. 11, 055063 (2009). 

41. V. Braude and E. B. Sonin, Phys. Rev. B 74, 064501 (2006). 

42. C. Azimonte, J.C. Cezar, E. Granado, Q. Huang, J.W. Lynn, J. CP. Campoy, J. 
Gopalakrishnan, and K. Ramesha, Phys. Rev. Lett. 98, 017204 (2007). 

43. S.W. Lovesey, K. S. Knight, and D. S. Sivia, Phys. Rev. B 65, 224402 (2002). 

44. I.V. Solovyev, PH. Dederichs, and I. Mertig, Phys. Rev. B 52, 13419 (1995). 

45. D. Xiao, Y. Yao, Z. Fang, and Q. Niu, Phys. Rev. Lett. 97, 026603 (2006). 

46. G. Autes, C. Barreteau, M. C. Desjonqueres, D. Spanjaard, and M. Viret, Europhys. 
Lett. 83, 17010 (2008). 

47. T. Koide, H. Miyauchi, J. Okamoto, T. Shidara, A. Fujimori, H. Fukutani, K. 
Amemiya, H. Takeshita, S. Yausa, T. Katayama, and Y. Suzuki, Phys. Rev. Lett. 
87, 257201 (2001). 

48. M.R. Geller and G. Vignale, Physica B 212, 283 (1995). 



Theory of Orbital Magnetization in Solids 29 



49. O. Entin-Wohlman, Y. Imry, A. G. Aronov, and Y. Levinson, Phys. Rev. B 51, 11584 
(1995). 

50. I.V. Solovyev, Phys. Rev. Lett. 95, 267205 (2005). 

51. J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1975). 

52. F. Mauri and S.G. Louie, Phys. Rev. Lett. 76, 4246 (1996). 

53. F. Mauri, B. G. Pfrommer, and S.G. Louie, Phys. Rev. Lett. 77, 5300 (1996). 

54. C.J. Pickard and F. Mauri, Phys. Rev. Lett. 88, 086403 (2002). 

55. D. Sebastiani and M. Parrinello, J. Phys. Chem. A 105, 1951 (2001). 

56. D. Sebastiani, G. Goward, I. Schnell, and M. Parrinello, Computer Phys. Commun. 
147, 707 (2002). 

57. L.L. Hirst, Rev. Mod. Phys. 69, 607 (1997). 

58. R. Resta, Rev. Mod. Phys. 66, 899 (1994). 

59. R. Wu, First Principles Determination of Magnetic Anisotropy and Magnetostriction 
in Transition Metal Alloys, Lecture Notes in Physics 580 (Springer, Berlin, 2001). 

60. S. Sharma, S. Pittalis, S. Kurth, S. Shallcross, J. K. Dewhurst, and E. K. U. Gross, 
Phys. Rev. B 76, 100401(R) (2007). 

61. F.D.M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

62. K. Ohgushi, S. Murakami, and N. Nagaosa, Phys. Rev. B 62, R6065 (2000). 

63. T. Jungwirth, Q. Niu, and A. H. MacDonald, Phys. Rev. Lett. 88, 207208 (2002). 

64. S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, 1348 (2003). 

65. Y. Yao, L. Kleinman, A. H. MacDonald, J. Sinova, T. Jungwirth, D.-S. Wang, E. 
Wang, and Q. Niu, Phys. Rev. Lett. 92, 037204 (2004). 

66. T. Thonhauser, D. Ceresoli, D. Vanderbilt, and R. Resta, Phys. Rev. Lett. 95, 137205 
(2005). 

67. N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 (1997). 

68. C. Brouder, G. Panati, M. Calandra, C. Mourougane, and N. Marzari, Phys. Rev. 
Lett. 98, 046402 (2007). 

69. L. He and D. Vanderbilt, Phys. Rev. Lett. 86, 5341 (2001). 

70. T. Thonhauser and D. Vanderbilt, Phys. Rev. B 74, 235111 (2006). 

71. R. Resta, D. Ceresoli, T. Thonhauser, and D. Vanderbilt, ChemPhysChem. 6, 1815 

(2005) . 

72. O. Gat and J.E. Avron, Phys. Rev. Lett. 91, 186801 (2003). 

73. D. Ceresoli, T. Thonhauser, D. Vanderbilt, and R. Resta, Phys. Rev. B. 74, 024408 

(2006) . 

74. I. Souza and D. Vanderbilt, Phys. Rev. B 77, 054438 (2008). 

75. D. Xiao, J. Shi, and Q. Niu, Phys. Rev. Lett. 95, 137204 (2005). 

76. X. Wang, J.R. Yates, I. Souza, and D. Vanderbilt, Phys. Rev. B 74, 195118 (2006). 

77. L. E. Reichl, A Modern Course in Statistical Physics (University of Texas Press, 
Austin, 1980). 

78. Q. Niu, Berry Phase and the Anomalous Hall Effect, talk presented at the 17 th 
Annual Workshop on Recent Developments in Electronic Structure Methods (Cornell 
University, June 24, 2005, Rhaca, New York). Available at |http :"7 /es05 .physics . 
Cornell . edu/schedule . shtml 

79. G. Sundaram and Q. Niu, Phys. Rev. B 59, 14915 (1999). 

80. J. Shi, G. Vignale, D. Xiao, and Q. Niu, Phys. Rev. Lett. 99, 197202 (2007). 

81. G. Vignale and M. Rasolt, Phys. Rev. B 37, 10685 (1988). 

82. N. Sai, K. M. Rabe, and D. Vanderbilt, Phys. Rev. B 66, 104108 (2002). 

83. I. Souza, J. Iniguez, and D. Vanderbilt, Phys. Rev. B 69, 085106 (2004). 

84. S. Baroni, S. Gironcoli, A. Dal Corso, and P. Giannozzi, Rev. Mod. Phys. 73, 515 
(2001). 



30 T. Thonhauser 



85. P.E. Blochl, Phys. Rev. B 50, 17953 (1994). 

86. See http://www.gipaw.net 

87. C.J. Pickard and F. Mauri, Phys. Rev. B 63, 245101 (2001). 

88. G. Schreckenbach and T. Ziegler, J. Phys. Chem. A 101, 3388 (1997). 

89. D. Ceresoli and R. Resta, Phys. Rev. B 76, 012405 (2007). 

90. R. Car and M. Parrinello, Phys. Rev. Lett. 55, 2471 (1985). 

91. R. Resta, in Berry Phase in Electronic Wavef unctions, Troisieme Cycle Lecture Notes 
(Ecole Polytechnique Federale de Lausanne, Switzerland, 1996). 

92. M. G. Lopez, D. Vanderbilt, T. Thonhauser, I. Souza, Phys. Rev. B, under review 
(2011). 

93. I. Souza, N. Marzari, and D. Vanderbilt, Phys. Rev. B 65, 035109 (2001). 

94. P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

95. R. M. Martin, Electronic Structure, Basic Theory and Practical Methods (Cambridge 
University Press, Cambridge 2004). 

96. M. Head-Gordon, Quantum Chemistry: Standard Methods and New Frontiers in 
Wavefunction Theory, talk presented at the KITP conference From Basic Con- 
cepts to Real Materials (University of California, November 5, 2009, Santa Bar- 
bara, California). Available at |http : / /online . kitp . ucsb . e du/online/excitcm_ 
c09/headgordon 

97. J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

98. J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

99. S.K. Ghosh and A. K. Dhara, Phys. Rev. A 38, 1149 (1988). 

100. G. Vignale, Phys. Rev. B 70, 201102(R) (2004). 

101. E. Runge and E.K. U. Gross, Phys. Rev. Lett. 52, 997 (1984). 

102. P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, 
G. L. Chiarotti, M. Cococcioni, I. Dabo, A. Dal Corso, S. de Gironcoli, S. Fabris, G. 
Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin- 
Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. Paulatto, 
C. Sbraccia, S. Scandolo, G. Sclauzero, A. P. Seitsonen, A. Smogunov, P. Umari, 
and R. M. Wentzcovitch, J. Phys.: Condens. Matt. 21, 395502 (2009). Available at 

: //www. quantum- espresso . org 

103. D. Bonnenberg, K. A. Hempel, H.P.J. Wijn, 1.2.1.2.4 Atomic Magnetic Moment, 
Magnetic Moment Density, g and g' Factor, in 1.2.1.2 Magnetic Properties, ed. H. P. 
J. Wijn. Part of Landolt-Bornstein - Group III Condensed Matter, Numerical Data 
and Functional Relationships in Science and Technology, Vol. 19a: 3d, 4d and 5d 
Elements, Alloys and Compounds. 

104. P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and J. Luitz, WIEN2k, An 
Augmented Plane Wave + Local Orbitals Program for Calculating Crystal Properties. 
Available at http : / /www . wien2k . at 

105. Y. Yao, First Principles Calculations of Orbital Magnetization - Preliminary Results, 
talk presented at the CECAM conference Orbital Magnetization in Condensed Matter 
(EPFL, June 15, 2009, Lausanne, Switzerland). Abstract available at http: //www. 
cecam. org/workshop-4-303 . html?presentation_id=3575 

106. 1. 1. Rabi, J. R. Zacharias, S. Millman, and P. Kusch, Phys. Rev. 53, 318 (1938). 

107. Encyclopedia of NMR, edited by D. M. Grant and R. K. Harris (Wiley, London, 1996). 

108. W. Kutzelnigg, U. Fleischer, and M. Schindler, NMR Basic Principles and Progress 
(Springer, Berlin, 1990). 

109. D. Ceresoli, N. Marzari, M. G. Lopez, and T. Thonhauser, Phys. Rev. B 81, 184424 
(2010). 

110. T. Thonhauser, D. Ceresoli, and N. Marzari, Int. J. Quant. Chem. 109, 3336 (2009). 



Theory of Orbital Magnetization in Solids 31 



111. G. Kresse, M. Marsman, and J. Furthmiiller, Vienna Ab Initio Simulation Package, 
VASP the Guide. Available at http://cms.mpi.univie.ac.at/vasp 

112. X. Gonze, B. Amadon, P.-M. Anglade, J.-M. Beuken, F. Bottin, P. Boulanger, F. 
Bruneval, D. Caliste, R. Caracas, M. Cote, T. Deutsch, L. Genovese, Ph. Ghosez, 
M. Giantomassi, S. Goedecker, D. R. Hamann, P. Hermet, F. Jollet, G. Jomard, 
S. Leroux, M. Mancini, S. Mazevet, M.J.T. Oliveira, G. Onida, Y. Pouillon, T. 
Rangel, G.-M. Rignanese, D. Sangalli, R. Shaltaf, M. Torrent, M. J. Verstraete, G. 
Zerah, and J. W. Zwanziger, Comp. Phys. Commun. 180, 2582 (2009). Available at 
http: //www. abinit . org 

113. G. te Velde, F.M. Bickelhaupt, S.J. A. van Gisbergen, C. Fonseca Guerra, E.J. 
Baerends, J. G. Snijders, and T. Ziegler, J. Comp. Chem. 22, 931 (2001). Available 
at http://www.scm.com 



